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Abstract. It has previously been shown heuristicahy that the topology of the 
Universe affects gravity, in the sense that a test particle near a massive object 
in a multiply connected universe is subject to a topologically induced acceleration 
that opposes the local attraction to the massive object. It is necessary to check if 
this effect occurs in a fully relativistic solution of the Einstein equations that has 
a multiply connected spatial section. A Schwarzschild-like exact solution that is 
multiply connected in one spatial direction is checked for analytical and numerical 
consistency with the heuristic result. The (slab space) heuristic result is found 
to be relativistically correct. For a fundamental domain size of L, a slow-moving, 
negligible-mass test particle lying at distance x along the axis from the object of mass 
M to its nearest multiple image, where GM/c^ <^ x L/2, has a residual acceleration 
away from the massive object of 4:({3)G{M/ L^) x, where C(3) is Apery's constant. For 
M ^ and L ^ 10 to 20h^^ Gpc, this linear expression is accurate to ±10% 

over 3h~^ Mpc<x<2h^^ Gpc. Thus, at least in a simple example of a multiply 
connected universe, the topological acceleration effect is not an artefact of Newtonian- 
like reasoning, and its linear derivation is accurate over about three orders of magnitude 
in X. 
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1. Introduction 



The topology of spatial sections of the U niverse has be e n of intere s t sin c e the 
foundations of r elativistic c osmol ogy 
Lemaitre. 



19271 : lR,obertsonl . \193$ 



le.g 
However 



de Sitter 



1917 



Friedmann 



1923 



1924 



astronomical measurements aimed at 



determining spatial topology have long suffered a theoretical disadvantage compared to 
curvature measurements. The latter are tightly related to another type of astronomical 
measurement (matter-energy density) via a well-established physical theory: general 
relativity. Some elementary work that might lead to a physical theory of cosmic 
topology has been carried out , by co mparing certain characteristics of different manifolds 



(ISeriul . Il996t lAnderson et al.l . l2004j ) and by explorations of so me elements of topology 



change in quantum gravity (e.g., iDowker and Suryal . Il998l ). but these remain very 
distant from astronomical observations. Thus, most of the empirical work of the 
past few decades has been carried out with no theoretical constraints, apart from 
assuming Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmological models. For 
recent empirical analyses of the case for and against a multiply connected positively 
curved spatial section, in partic ular the Poincare dodecahedral space S^/J*, see 
Roukema and KazimierczakI (l201ll ) and references therein. Theoretical developments 
could offer the possibility of new astronomical tests for cosmic topology. 

For perfectly homogeneous solutions of the Einstein equations, i.e. FLRW models, 
if the sign of the curvature is determined empirically, only the covering space H^, M^, or 
S^, i.e. an apparent space containing many copies of the fundamental domain, can be 
inferred. The 3-manifold of comoving space itself, i.e. H^/F, R^/F, or S^/F, respectively, 
for a fundamental group (holonomy group) F, is not implied by the curvature. However, 
heuristic, Newtonian-like calculations have recently shown that the presence of a single 
inhomogeneity is enough to, in pr i nciple, distinguish different possible s paces of the same 
fixed curvature ( iRoukema et al.l . 120071 : iRoukema and Rozahski l2009l ). The effect is a 
long-distance, globally induced acceleration term that opposes the local attraction to 
a massive object. Curiously, the effect, (hereafter, "topological acceleration''^) is much 
weaker in the Poincare dodecahedral space than in other spaces, hinting at a theoretical 
reason for selecting this space independent of the observational reasons published six 
years earlier. Other calculations of globally induced effec ts of the spatial topology of a 
fixed FLRW background include llnfeld and Schildl (119461 ) and iBernui et al.l (119981 ) . 

Astronomical tests for cosmic topology are almost always based on the fact that 
photons from a p oint in comoving space can arrive at the observer by multiple paths (e.g. 
Roukemal . l2002l ) . In contrast, topological acceleration has the potential to provide tests 
of galaxy kinematics that could provide experimental constraints on cosmic topology 
independent of tests based on multiple imaging (of extragalactic objects or of cosmic 
microwave background temperature fluctuations). However, the calculations of the 

I See Lemaitr j (1931) for an incomplete translation of iLemaitre (1927). 

§ The implied meaning is "acceleration induced globally by cosmic topology" ; a term used in earlier 
work is "residual gravity" . 
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effect (iRoukema et al.l . 120071 : iRoukema and Rozanskil . |2009| ) did not use full solutions 
of the Einstein equations. The FLRW metric is an exact solution of the Einstein 
equations, but it can only by applied to the real Universe as a heuristic guide, since 
galaxies (inhomogeneities) exist. When this heuristic guide is used to interpret recent 
observational constraints, an accelerating scale factor and dark ener gy are inferred . That 
is, the latter appear to be artefacts of a heuristic appro ach (e.g. 



Buchert and Carfora . 



2003 



Celerier et al 



2010 



Wiegand and Buchertl . l2010l ). Similarly, since a heuristic. 



Newtonian-like approach was used to calculate the topological acceleration effect, the 
latter could, in principle, also be an artefact. Thus, it is important to check whether 
the effect really exists in the relativistic case of a (3 + l)-dimensional spacetime 
with multiply connected spatial sections. Indeed, in reviews of cosraic topology (e.g 



EUisl. 



1971 



Lachieze-Rey and Luminetl . Il995l : iLuminet and Roukemal . 11999 



Levinl . 12002 



Reboucas and Gomerol . 20041 ) . the use of the exact (locally homogeneous) FLRW models 



has generally led to the inference that in a classical general relativistic model, spatially 



methods, see e.g.. 


Uzan et al. 


1999: 


Roukema and Luminet. 


1999: 


Bento et al.. 


2006: 


Reboucas and Alcanizl. 


2006 


; 1 


3ento et al.[ 


2006: 


Reboucas et al.. 


2006 


). This inference 



only becomes incorrect when the universe model contains at least one inhomogeneity. 

Strictly speaking, from the point of view of differential geometry alone, astronomical 
tests for cosmic topology based on multiple imaging are degenerate. They do not 
distinguish a multiply connected space from a simply connected universal cover that 
happens to be populated by physically distinct objects located at appropriate positions 
in a perfectly regular tilingljjj Similarly, what we refer to here as "topological 
acceleration" could be interpreted as an effect of a perfectly regular tiling in a simply 
connected covering space. In this paper, we adopt the Schwarzschildian approach to 
lifting the degeneracy, 'We would be much happier with the view that these repetitions 
are illusory, that in reality space has peculiar connection properties so that if we leave 
any one cube through a side, then we immediately reenter it through the opposite side," 
(transl., ISchwarzschildl . Il900l : IStewart et al.l . Il998[ l This is similar to the interpretation 
of what are normally claimed to be multiple, "gravitationally lensed" images of a single 
physical object as representing a single object rather than multiple objects that happen 



tomimic what is expected from a gravitational lensing hypothesis (e.g., lAdam et al. 
1989h . 



Here, we compare the heuristic calculation of IRoukema et al.l ( 120071 ) for T^ (slab 
space, i.e. S^ x R^) to a limiting case of an exact, Schwarzschild-like solution 
which has a T^ spatial section outside of the event horizon. This is a simpler 3- 
mar iifold than those that seern to b e favoured by observations (as stated above. 



see 



Roukema and Kazimierczakl . 1201 ll . and references therein), with only one short 



dimension, and it only contains one massive object. Checking the existence of topological 

II If a spatially closed, timelike path (world-line) were allowed around a non-contractible spatial loop 
in a given universe model, then a physical voyage by an observer around the loop would provide an 
experimental method to overcome this degeneracy. 
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Figure 1. Event horizon of massive object of mass M (dark line segment) and 
test particle (dot) in slab space of fundamental domain length L. In a 3-section 
of constant coordinate time t, the event horizon of the massive object is spherically 
symmetric in terms of the metric, but in the Weyl coordinates used here, it is a line 
segment extending from — M to +M in the x direction, where G ^ c — 1, and is of 
zero size in the radial coordinate p. The test particle is displaced a small distance from 
the massive object, i.e. x <C L/2, but well outside the event horizon, i.e. < M <C a; 
(zoom) [Eq.©]. 



acceleration in this relativistic case does not guarantee corresponding results in more 
complicated cases, but it does open the prospect of generalisations. In Sect. 12.11 we 
list the simplifying assumptions for considering a test particle distant from the event 
horizon but much closer to the "first" copy of a massive object than to an "adjacent" 
copy. The metric in Weyl coordinates and related functions (IKorotkin and Nicolail . 
I994J ) are given in Sect. 12.21 for an analytical derivation and the numerical method is 
stated in Sect. 13.21 The analytical and numerical calculations are presented in Sect. 13.11 
and Sect. 13. 2[ respectively, and a conclusion is given in Sect. HI 



2. Method 



For the Newtonian-like calculation. Fig. 3 of iRoukema et al.l (120071 ) illustrates the 
situation of a test particle at distance x from a point object of mass M in T^, lying along 
the spatial geodesic of length L from the massive o bject to itself. Definin g e := x/L 
gives the topologically induced acceleration (Eq. (7). iRoukema et al.l . 120071 ) 



X 



when e ^ 1, adopting G 



{jL - x)2 
= 1. 



ijL + x)^ 



4.8 



Me 
1? 



(1) 



2.1. Assumptions 



A Schwarzschild-like solution of the Einstein equations has an event horizon, but 
behaviour close to the event horizon (or inside it) is not of interest here. Also, if the 



Relativistic topological acceleration 



5 



topological acceleration effect were to be analysed for its effects on galaxy dynamics, 
peculiar velocities would typically be of the order of 10~^c and very rarely exceed 10~^c. 
As in the earlier derivation of Eq. ([T]) , we consider a test particle lying along the closed 
spatial geodesic joining the massive object to itself, i.e. having Weyl radial coordinate 
p = 0. 

Thus, in addition to adopting the conventions G = c = 1, a (—,+,+,+) metric, 
Greek spacetime indices 0,1,2,3 and Roman space indices 1,2,3, our assumptions 
become: 

< M < a; < L/2, (2) 
dx* ^ dx* ^ ^ dt 

_«1 ^ ^«l<^. (3) 
p = 0, (4) 

where r is the proper time along the test particle's world line x°(r), and t = x^. The 
assumption of a low coordinate velocity implies a low 4-velocity spatial component 
[Eq. ©]. 

2.2. Analytical approach 



Korotkin and Nicolail (jl994j ) presented a family of multiply connected exact solutions 
of the Einstein equations that includes a Schwarzschild-like solution with T^ spatial 
sections outside of the event horizon. Figure [1] shows this solution in Weyl coordinates 
X and p. Using Weyl coordinates and the Ernst potential, the metric is 

ds' = -e'dt^ + e-" [e^'idx^ + dp^) + pM^^] (5) 



from Eq. (9) of lKorotkin and Nicolai where k and u relate to an Ernst potential 



e defined on ^ := x + ip, with 



v/(x - M)2 + p2- + ^{x + M)2 +"7 - 2M 
v/(x - M)2 + p2 + ^{x + MY + p2 + 2M ^' 
2M 

loq := Ineo, oq := 0, a^yo := -^-y (7j 

oo 

^ix,p)-= ^ Mx + jL,p) + aj], e:=e'^ (8) 

j=-oo 

a,..2,p|£|l. (9) 



from Eqs (13), (12), (5) and (7) of lKorotkin and Nicolail JiQqJ I+I The Ernst potential 



is real for this solution, so Eq. (|9]) gives 

d,k = iip^ = lipm^nef = hpid^uY. (10) 



% In the online version 1 of arXiv : gr-qc/9403029, the time component has an obvious typological 
error in the sign. 

+ Positive square roots are implied in Eq. 
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Thus, at p = 0, we have d^k = 0, i.e. /c is a constant along p = 0. We choose 
k{p = 0) = 0, since any other value just implies a rescaling of units. Thus, the metric 
on p = is 

ds^ = -e'dt' + e-" {dx^ + dp2 + pM02) . (11) 

Defining the tangent vector V by components V"' := dx"/dr, the geodesic equation 
for the test particle is Vy V = 0, i.e. 

d^x" dx^dxT_ 

17^ + ^ ^^17d7""' 
where F"^^ are the Christoffel symbols of the second kind|3 Since a (— , +, +, +) 
convention is used and low velocities are assumed [Eq. ([3])], coordinate and proper 
time are related by 

dT' = \goo\{dxy = e-dt', (13) 

giving the coordinate acceleration 

^ = ^^ = P— (U) 
dt2 dt2 dr2' ^ ' 

Given the conditions in Eqs. ([2]), (jS]), and (jl]) and the Ernst potential expressions in 

Eqs. f lTTj) . ([6]), ([7]), and ([8]), conversion to metric units of x and t would only modify this 

to second order. In Sect. 13. H d^x/dt^ is evaluated. 

2.3. Numerical approach 

The Ernst potential expressions [Eqs ([6]), ([7]), and ([8])] are used for numerical evaluation 
of Eq. ( 120|) (below) by finite differencing over small intervals. A typical cluster of 
galaxies has a mass of ~ IO^^Mq, i.e. M = 4.78 pc in length units, and the most 
massive clusters have masses up to about 10^^Mr;,fl Observational estimates of L are 



in the range 10 to 20h ^ Gpc (e.g. iRoukema and Kazimierczakl . 1201 ll . and references 



therein). Thus, typical scales of M/L of interest should be M/L ~ 10~ to 10^ . Since 
M/x,M/L, and x/L are small, the use of 53-bit-significand double-precision floating- 
point numbers is replaced by arbitrary precision arithmetic, for 10"^° < M/L < 10"®. 
These calculations are presented in Sect. 13.21 

3. Results 

3.1. Coordinate acceleration 

Given that the test particle is far from the "close" copy of the event horizon and from 
"distant" copies [Eq. (|2])], we have |r"^^| ^ 1. Together with Eq. (|3]), this implies that 

* Symmetric definition. 

ji In the Weyl coordinate x, the event horizon occurs along |a;| < M, p — (Fig. [1] which may 
seem surprising g iven t hat th e event horizon in Schwarzschild coordinates occurs at r = 2M. See 



Frolov and Frolovl (e.g., l2003l Sect. II. B. and references therein) for the relation between Weyl and 



Schwarzschild coordinates. 
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most of the terms in the sum in Eq. f ll2p are small terms of second order or higher, 
leaving 



+ r 



00 



dx° 

d7 



0. 



Since the metric is static (and diagonal), we have 

r"oo = ^^"^(^o^Ao + <9o^oA - <9a^oo) = -^^"^<9a^oo- 
Using Eq. (fT3|) . Eq. ( ITSj) becomes 

d^x" _ 1 ^ 
dr^ 2 1 5(00 1 

Since the metric is diagonal, we have g^^ = (gxx)^^ = e^'^ and 
d^x^ 1 



I.e., 



dr2 2 1 500 
d^x 1 g 



■digoo, 



dr2 2 |5it| 
Thus, from Eq. f lT4|) . 



1 e'^ 

9.5« = 2^^-^-^' 



2 



d^x U 

-rrr = <^x^- 

dt2 2 



(15) 



(16) 



(17) 



:i8) 



(19) 



(20) 



Subst ituting Eq ([7]) in Eq. (|8^) [corresponding to (12) and (5) of iKorotkin and Nicolai 
(I1994J ). respectively] 

2M' 



u{x,p) = wo{x,p) + ^ 

j=-oojyo 



uoix + jL,p) + 



L\j\ 



(21) 



This is the convergent solution found by iKorotkin and Nicolail ( 19941 ). Dropping the p 
dependence since we are interested in the p = axis, the derivative can be written using 
Eq. (CD 



dxUj{x) = dx 



lneo{x)+ ^ ln£:o(a; + j-^^) 

j=-oojVO 



^0:^0(3^) ^ dx€o{x + jL) ^ y> dxeo{x-jL) 
dxSoix + jL) ^ dxSoijL - x) 



dxSoix 



^o(a;) jr[ £o{x+jL) ^ soijL-x) 



1= 

00 



fa dxeo{x) + ^ dxeoiJL + a;) - ^ dxSoUL - x), 

(22) 

where the third equality follows from using Ox'Sq^—x') = —dx'eo{x') [cf. Eq. ([6])], and 
the approximation follows from using Eqs ([2]), (jlj), and j > 1 in Eq. 
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1e-05 0.0001 0.001 0.01 0.1 1 

X 



Figure 2. Topological acceleration di^x/dt^ + M/x^ as a function of coordinate 
distance x from the "first" copy of the massive object of mass M , along the spatial 
geodesic joining the massive object to itself, using Eqs ([B]), ©, and ([5]) in Eq. (pOj) . 
for L = 1. The curves show masses from M = 10""'^° to 10"^, from bottom to top, 
respectively, as labelled. The sloping lines show 4C(3)(Af/i^)a;. The vertical line at 
X = 0.5 shows the halfway point to the "next" copy of the massive object. The vertical 
lines down to and up to the M = 10"^'' curve show the range over which the linear 
approximation is valid to within ±10%. 

Table 1. Minimum and maximum distance x for which the 4(^(3) x 
approximation is accurate to within 10% (cf. Fig. [2]).'^ 



M 


10-10 


3 X 10-1° 


10-9 


3 X 10-9 


10-8 




-3.94 


-3.77 


-3.54 


-3.37 


-3.14 




-0.63 


-0.63 


-0.63 


-0.63 


-0.63 


A logio 


3.32 


3.14 


2.91 


2.74 


2.51 



^For convenience, L = 1. 
^Range in log^o x. 



For < M < x', Eq. (jSD simplifies to 
x' -M 



eo{x',0) 



x' + M 

Again defining e := x/L, we then have 

dxOj{x) 

2M ^ 2M 



(23) 



2M 



X + M)2 (jL + X + M)2 ^ {jL -x + M) 



2M 2M 

+ 



x^ 



L2 



2M 2M 



oo 



^ 1 / 2e\ ^ 1 f 2e 
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2M 



8M ^ 1 



^2M_SCi0_^ (24) 

where Eq. ([2]) is applied in the second and third hnes, and ({3) is the Riemann zeta 
function evaluated at an argument of 3, i.e. Apery's constant. 
Hence, using Eq. (12UI) . 



d^x ^ M ^ 4C(3)M X 



M 4.8M X 



dt^ ' x2 ' L2 L " x2 ' L2 L 
i.e. the coordinate acceleration is the Newtonian acceleration towards the local copy of 
the massive object plus an oppo sing topological a cceler ation term matching Eq. ([T]), i.e. 
the expression derived earlier in iRoukema et al.l ( 120071 ) . 



3.2. Numerical check 

For 10"^° < M < 10"^, L = 1, \j\ < 14, and linear offsets of Sx = ±10"^^ around each 
X value, numerically convergent results were found using 100 bits in the significand. 
Figure [2] and Table [T] show that in this case, the linear approximation for the residual 
acceleration, 4^(3) (M/L'^) x, is accurate to within 10% over about three orders of 
magnitude. Near the halfway point and at x > 0.5, i.e. closer to the "next" image 
of the massive object than the "local" copy, the residual term defined in relation to 
the acceleration towards the "local" copy is clearly no longer small — it diverges as the 
test particle approaches the "next" copy. Thus, for L ~ 10 to 20h~^ Gpc, the linear 
expression should be accurate on scales of 3h~^ Mpc ^ x <, 2h^^ Gpc in a slab-space 
universe containing just one cluster of about 



4. Conclusion 

The T^ (slab space) heuristic result is found to exist in the limit of the Schwarzschild- 



l ike, e xact, slab-space solution of the Einstein equations found by lKorotkin and Nicolai 



( 1l994j ). with the same linear topological acceleration term as for the simpler calculation. 
For a fundamental domain size of L, a slow- moving [Eq. ([3])], low-mass test particle lying 
at distance x ^ L/2 along the axis [p = 0, Eq. (jl])] from the object of positive mass 
GM/c^ ^ X [Eq. (|2])] to its nearest multiple image, the additional acceleration away 
from the massive object is A({3)G{M/ L^) x [Eq. ( |25|) ]. For a topological acceleration 
term away from a cluster of ~ 10^^ Mq and a Universe of siz 60 L ~ 10 to 20/1-1 Gpc, 
the linear expression should be accurate over 3h~^ Mpc ^x ^ 2h~^ Gpc (Fig. |2l column 
with M = 3 X 10-1° in Table [I]). Thus, at least in a simple example of a multiply 

ft Mo re formally, twice the injectivity radius, or twice the in-radius (Fig. 10, iLuminet and Roukemal . 



1999L and references therein). 
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connected universe, the topological acceleration effect is not an artefact of Newtonian- 
like reasoning. The linear term exists as an analytical limit of the relativistic model and 
as a good numerical approximation to it over several orders of magnitude of x. 

This qualitatively suggests that the Newtonian-hk e derivations of the topologica l 



acceleration effect for well-proportioned FLRW models (iRoukema and Rozahskil . l2009l ). 
in which the Poincare dodecahedral space is found to be uniquely well-balanced, might 
also be valid limits of fully relativistic Schwarzschild-like solutions with these spatial 
sections. Other extensions of the present work within the slab-space solution could 
include consideration of test particles with high coordinate velocities [violating Eq. ([3])], 
or expanding (FLRW-like) solutions. 

Although we have shown that topological acceleration exists theoretically, the 
effect is clearly weak over the range for which the linear expression is accurate. It 
is not clear how easy it might be to test the effect observationally. A possible method 
might follow from considering the effect as a feedback effect from global anisotropy to 
local anisotropy. Multiply connected FLRW models have usually been thought to be 
locally isotropic (the metric on a spatial section is isotropic) and globally anisotropic 
(some global tests differ depending on the chosen spatial direction). Here, we have 
confirmed that global anisotropy can imply a local acceleration effect. The effect is 
clearly anisotropic. The integration of this weak effect over a long period of time for 
individual particle trajectories or collapsing density perturbations might, in principle, 
give anisotropic long-term statistical behaviour that is strong enough to be detectable. 
Thus, the global anisotropy of a multiply connected space implies (at least in the case 
considered here) a weak, locally anisotropic dynamical effect. Whether or not the time- 
integrated effect could be strong enough to detect in observations of galaxy kinematics 
or cosmic web structure is a question worth exploring in future work. 
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